IN THE CLAIMS : 

Please amend the claims as follows: 

1 . (Currently Amended) A method of adaptive direct volume rendering, 
comprising 

fragmenting a sampled 3-D dataset of a scalar field into a plurality of sub- 
volumes of different sizes, each sub-volume associated with a set of data value parameters 
characterizing the data value distribution of the scalar field within the sub-volume; 

defining an opacity transfer function that is dependent upon data values of the 
scalar field and an illumination model; 

selectively casting a plurality of rays from a 2-D image plane towards the 
sampled dataset, each ray having an initial ray energy and a cross-section; 

for each ray cast from a selected location on the 2-D image plane, 

electing a subset of the plurality of sub- volumes for interacting with 

the ray; 

estimating the ray energy reflected by each sub-volume of the subset 
using the opacity transfer function and the illumination model; and 

summing the reflected ray energy as a pixel value at the selected 
location on the 2-D image plane; and 

estimating pixel values at other locations on the 2-D image plane using the 
pixel values at the selected locations^ 

wherein the step of selectively casting comprises; 
subdividing the 2-D image plane into a plurality of sub-planes; 
for each of the plurality of sub-planes, 

casting four rays from the four corners of the sub-plane and estimating 
a pixel value at each corner; 

calculating a maximum pixel value variation within the sub-plane; and 
recursively subdividing the sub-plane into multiple child sub-planes of 
smaller sizes by casting a ray from the center of the sub-plane until the maximum pixel value 
variation of the sub-plane is below a predefined imaging error threshold , 

2. (Original) The method of claim 1 , wherein 

the step of fragmenting a sampled 3-D dataset includes 

fragmenting the 3-D dataset into eight sub-volumes; and 
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for each sub-volume, recursively fragmenting it into eight smaller sub- 
volumes until the size of the smallest sub-volumes reaches a predefined size limit. 

3. (Original) The method of claim 2, wherein the predefined size limit is a sub- 
volume comprising 2x2x2 3-D cells and the eight corners of each cell are associated with 
eight data values of the scalar field. 

4. (Original) The method of claim 3, wherein the data value at a location within 
the cell is tri-linearly interpolated using the eight data values at the eight corners of the cell. 

5. (Original) The method of claim 1, wherein the set of parameters include a 
maximum, an average value, and a minimum data value of the scalar field within the sub- 
volume. 

6. (Original) The method of claim 1, further comprises 

constructing an octree comprising a root node, a plurality of intermediate 
nodes, and a plurality of leaf nodes; 

associating the root node with the 3-D dataset; 

associating each of the plurality of leaf nodes with a smallest sub-volume from 
the plurality of sub-volumes; and 

associating each of the plurality of intermediate nodes with a sub-volume from 
the plurality of sub-volumes that is larger than the smallest sub-volume. 

7. (Canceled) 

8. (Currently amended) The method of claim [[7]] i, wherein the maximum pixel 
value variation within the sub-plane is defined as the maximum deviation of pixel values at 
the four corners of the sub-plane from an average pixel value of the sub-plane. 

9. (Currently amended) The method of claim [[7]] i, wherein the predefined 
imaging error threshold is modulated by an image rendering speed provided by a user, a 
distance to an edge of an object embedded in the 3-D dataset, and a difference between a 
pixel value estimated from an adaptive ray casting and a pixel value estimated from a bi- 
linear interpolation. 

10. (Currently amended) Tho method of claim 1, A method of adaptive direct 
volume rendering, comprising 

fragmenting a sampled 3-D dataset of a scalar field into a plurality of sub- 
volumes of different sizes, each sub- volume associated with a set of data value parameters 
characterizing the data value distribution of the scalar field within the sub-volume; 
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defining an opacity transfer function that is dependent upon data values of the 
scalar field and an illumination model; 

selectively casting a plurality of rays from a 2-D image plane towards the 
sampled dataset each ray having an initial ray energy and a cross-section; 

for each ray cast from a selected location on the 2-D image plane, 

selecting a subset of the plurality of sub-volumes for interacting with 

the ray; 

estimating the ray energy reflected by each sub-volume of the subset 
using the opacity transfer function and the illumination model; and 

summing the reflected ray energy as a pixel value at the selected 
location on the 2-D image plane; and 

estimating pixel values at other locations on the 2-D image plane using the 
pixel values at the selected locations 

wherein the step of selecting a subset of the plurality of sub-volumes for 
interacting with the ray includes 

identifying a largest sub-volume along the ray path and its corresponding 
maximum and minimum data values; 

checking if the opacity transfer function varies monotonically between the 
maximum and minimum scalar field values; 

if the function does not vary monotonically, recursively 

identifying a smaller sub- volume along the ray path and its 
corresponding maximum and minimum data values; and 

checking if the opacity transfer function varies monotonically between 
the maximum and minimum scalar field values of the smaller sub-volume; and 

if the function does vary monotonically, calculating the amount of ray energy 
reflected by the sub-volume during its interaction with the ray. 

1 1 . (Original) The method of claim 10, wherein two lookup tables are constructed 
for the opacity transfer function such that a forward lookup table contains the data value 
difference to a nearest local extreme of the opacity transfer function along the data value 
increasing direction and a backward lookup table contains the data value difference to a 
nearest local extreme of the opacity transfer function along the data value decreasing 
direction. 
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12. (Original) The method of claim 1 1 , wherein if the maximum data value of the 
sub-volume is smaller than the summation of the minimum data value of the sub-volume and 
its corresponding data value difference stored in the forward lookup table or the minimum 
data value of the sub- volume is larger than the difference between the maximum data value 
of the sub-volume and its corresponding data value difference stored in the backward lookup 
table, the opacity transfer function varies monotonically between the minimum and maximum 
data values. 

1 3 . (Currently amended) Tho method of claim 1, A method of adaptive direct 
volume rendering, comprising 

fragmenting a sampled 3-D dataset of a scalar field into a plurality of sub- 
volumes of different sizes, each sub- volume associated with a set of data value parameters 
characterizing the data value distribution of the scalar field within the sub- volume; 

defining an opacity transfer function that is dependent upon data values of the 
scalar field and an illumination model; 

selectively casting a plurality of rays from a 2-D image plane towards the 
sampled dataset, each ray having an initial rav energy and a cross-section; 

for each ray cast from a selected location on the 2-D image plane, 

selecting a subset of the plurality of sub-volumes for interacting with 

the ray; 

estimating the ray energy reflected by each sub-volume of the subset 
using the opacity transfer function and the illumination model; and 

summing the reflected ray energy as a pixel value at the selected 
location on the 2-D image plane; and 

estimating pixel values at other locations on the 2-D image plane using the 
pixel values at the selected locations 

wherein the step of estimating the ray energy reflected by each sub-volume of 
the subset includes 

estimating a maximum energy differential of the sub-volume; 

comparing the maximum energy differential against a predefined energy error 

threshold; 

if the maximum energy differential is above the predefined energy error 
threshold, recursively 
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selecting a smaller sub-volume along the ray path; and 
estimating a new maximum energy differential of the smaller sub- 
volume; and 

if the maximum energy differential is below the predefined energy error 
threshold, calculating the amount of ray energy reflected by the sub-volume using the 
illumination model. 

14. (Original) The method of claim 13, wherein the maximum energy differential 
depends on the opacity transfer function and the maximum, average, and minimum data 
values of the sub-volume. 

15. (Original) The method of claim 13, wherein the amount of ray energy 
reflected by the sub-volume depends on the length of ray path within the sub-volume, the 
opacity transfer function within the sub-volume, the average scalar field value of the sub- 
volume, and the local gradient vector of scalar field within the sub-volume. 

16. (Original) The method of claim 13, wherein if the sub-volume is a smallest 
sub-volume comprising 2x2x2 3-D cells, the smaller sub-volume is a 3-D cell within the 
smallest sub-volume, and if the opacity function does not vary monotonically within the cell, 
the 3-D cell is further divided into multiple sub-cells until the dimension of a smallest sub- 
cell reaches the cross-section of the ray. 

17. (Original) The method of claim 13, wherein if the sub-volume is a smallest 
sub-volume comprising 2x2x2 3-D cells, the smaller sub-volume is a 3-D cell within the 
smallest sub-volume, and if the opacity function varies monotonically within the cell, the 
maximum energy differential of the 3-D cell is estimated by dividing the maximum energy 
differential of the sub-volume by 2. 

18. (Original) The method of claim 13, wherein if the sub- volume is a smallest 
sub-volume comprising 2x2x2 3-D cells, the smaller sub-volume is a 3-D cell within the 
smallest sub- volume, and if the opacity transfer function varies monotonically within the cell 
and an iso-surface exists in the 3-D cell, the maximum energy differential of the 3-D cell is 
calculated using the eight data values at the corners of the 3-D cell and the opacity transfer 
function. 

19. (Original) The method of claim 13, wherein the predefined energy error 
threshold is modulated by an image rendering speed specified by a user and a zoom factor in 
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the case of parallel projection or a perspective angle and a perspective distance between the 
image plane and the 3-D dataset in the case of perspective projection. 

20. (Original) The method of claim 1, wherein 

the step of estimating pixel values at other locations on the 2-D image plane 

includes 

for each location, 

selecting four pixel values associated with four ray origins surrounding 

the location; and 

bi-linearly interpolating a pixel value at the location using the four 

pixel values. 

2 1 . (Currently Amended) An adaptive direct volume rendering system, 
comprising: 

one or more central processing units for executing programs; 
a user interface for receiving a plurality of volume rendering parameters; and 
an adaptive volume rendering engine module executable by the one or more central 
processing units, the module comprising: 

instructions for fragmenting a sampled 3-D dataset of a scalar field into a 
plurality of sub-volumes of different sizes, each sub-volume associated with a set of data 
value parameters characterizing the data value distribution of the scalar field within the sub- 
volume; 

instructions for defining an opacity transfer function that is dependent upon 
data values of the scalar field and an illumination model; 

instructions for selectively casting a plurality of rays from a 2-D image plane 
towards the sampled dataset, each ray having an initial ray energy and a cross-section; 

for each ray launched from a selected location on the 2-D image plane, 

instructions for selecting a subset of the plurality of sub-volumes for 
interacting with the ray; 

instructions for estimating the ray energy reflected by each sub- volume 
of the subset using the opacity transfer function and the illumination model; and 

instructions for summing the reflected ray energy as a pixel value at 
the selected location on the 2-D image plane; and 



1 -NY/1 998980.1 



-7- 



instructions for estimating pixel values at other locations on the 2-D 
image plane using the pixel values at the selected locations 

wherein the instructions for selectively casting a plurality of rays from a 2-D 
image plane include 

subdividing the 2-D image plane into a plurality of sub-planes; 
for each of the plurality of sub-planes, 

casting four rays from the four corners of the sub-plane and estimating 
a pixel value at each corner; 

calculating a maximum pixel value variation within the sub-plane; and 
recursively subdividing the sub-plane into multiple child sub-planes of 
smaller sizes by casting a ray from the center of the sub-plane until the maximum pixel value 
variation of the sub-plane is below a predefined imaging error threshold , 

22 . (Original) The system of claim 2 1 , wherein 

the step of fragmenting a sampled 3-D dataset includes 

fragmenting the 3-D dataset into eight sub-volumes; and 
for each sub-volume, recursively fragmenting it into eight smaller sub- 
volumes until the size of the smallest sub- volumes reaches a predefined size limit. 

23. (Original) The system of claim 22, wherein the predefined size limit is a sub- 
volume comprising 2x2x2 3-D cells and the eight corners of each cell are associated with 
eight data values of the scalar field. 

24. (Original) The system of claim 23, wherein the data value at any location 
within the cell is tri-linearly interpolated using the eight data values at the eight corners of the 
cell. 

25. (Original) The system of claim 21, wherein the set of parameters include a 
maximum, an average value, and a minimum data value of the scalar field within the sub- 
volume. 

26. (Original) The system of claim 2 1 , further comprises 
instructions for constructing an octree comprising a root node, a plurality of 

intermediate nodes, and a plurality of leaf nodes; 

instructions for associating the root node with the 3-D dataset; 

instructions for associating each of the plurality of leaf nodes with a smallest sub- 
volume from the plurality of sub-volumes; and 
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instructions for associating each of the plurality of intermediate nodes with a sub- 
volume from the plurality of sub-volumes that is larger than the smallest sub-volume. 

27. (Canceled) 

28. (Currently amended) The system of claim [[27]] 2j_, wherein the maximum 
pixel value variation within the sub-plane is defined as the maximum deviation of pixel 
values at the four corners of the sub-plane from an average pixel value of the sub-plane. 

29. (Currently amended) The system of claim [[27]] 21, wherein the predefined 
imaging error threshold is modulated by an image rendering speed provided by a user, a 
distance to an edge of an object embedded in the 3-D dataset, and a difference between a 
pixel value estimated from an adaptive ray casting and a pixel value estimated from a bi- 
linear interpolation. 

30. (Currently amended) Th e syst e m of claim 21, An adaptive direct volume 
rendering system, comprising: 

one or more central processing units for executing programs; 
a user interface for receiving a plurality of volume rendering parameters; and 
an adaptive volume rendering engine module executable by the one or more 
central processing units, the module comprising: 

instructions for fragmenting a sampled 3-D dataset of a scalar field into 
a plurality of sub-volumes of different sizes, each sub-volume associated with a set of data 
value parameters characterizing the data value distribution of the scalar field within the sub- 
volume; 

instructions for defining an opacity transfer function that is dependent 
upon data values of the scalar field and an illumination model; 

instructions for selectively casting a plurality of rays from a 2-D image 
plane towards the sampled dataset, each ray having an initial rav energy and a cross-section; 

for each ray launched from a selected location on the 2-D image plane, 
instructions for selecting a subset of the plurality of sub- 
volumes for interacting with the ray; 

instructions for estimating the rav energy reflected bv each sub- 
volume of the subset using the opacity transfer function and the illumination model; and 

instructions for summing the reflected ray energy as a pixel 
value at the selected location on the 2-D image plane; and 



1 -NY/1 998980.1 



-9- 



instructions for estimating pixel values at other locations on the 2-D 
image plane using the pixel values at the selected locations 

wherein the instructions for selecting a subset of the plurality of sub-volumes for 
interacting with the ray include 

identifying a largest sub-volume along the ray path and its corresponding 
maximum and minimum data values; 

checking if the opacity transfer function varies monotonically between the 
maximum and minimum scalar field values; 

if the function does not vary monotonically, recursively 

identifying a smaller sub- volume along the ray path and its 
corresponding maximum and minimum data values; and 

checking if the opacity transfer function varies monotonically between 
the maximum and minimum scalar field values of the smaller sub- volume; and 

if the function does vary monotonically, calculating the amount of ray energy 
reflected by the sub-volume during its interaction with the ray. 

3 1 . (Original) The system of claim 30, wherein two lookup tables are constructed 
for the opacity transfer function such that a forward lookup table contains the data value 
difference to a nearest local extreme of the opacity transfer function along the data value 
increasing direction and a backward lookup table contains the data value difference to a 
nearest local extreme of the opacity transfer function along the data value decreasing 
direction. 

32. (Original) The system of claim 3 1 , wherein if the maximum data value of the 
sub-volume is smaller than the summation of the minimum data value of the sub-volume and 
its corresponding data value difference stored in the forward lookup table or the minimum 
data value of the sub-volume is larger than the difference between the maximum data value 
of the sub-volume and its corresponding data value difference stored in the backward lookup 
table, the opacity transfer function varies monotonically between the minimum and maximum 
data values. 

33 . (Currently amended) Tho oyotom of claim 21 , An adaptive direct volume 
rendering system, comprising: 

one or more central processing units for executing programs; 

a user interface for receiving a plurality of volume rendering parameters; and 
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an adaptive volume rendering engine module executable by the one or more 
central processing units, the module comprising: 

instructions for fragmenting a sampled 3-D dataset of a scalar field into 
a plurality of sub-volumes of different sizes, each sub-volume associated with a set of data 
value parameters characterizing the data value distribution of the scalar field within the sub- 
volume; 

instructions for defining an opacity transfer function that is dependent 
upon data values of the scalar field and an illumination model; 

instructions for selectively casting a plurality of rays from a 2-D image 
plane towards the sampled dataset, each ray having an initial ray energy and a cross-section; 

for each ray launched from a selected location on the 2-D image plane, 
instructions for selecting a subset of the plurality of sub- 
volumes for interacting with the ray; 

instructions for estimating the ray energy reflected by each sub- 
volume of the subset using the opacity transfer function and the illumination model; and 

instructions for summing the reflected ray energy as a pixel 
value at the selected location on the 2-D image plane; and 

instructions for estimating pixel values at other locations on the 2-D 
image plane using the pixel values at the selected locations 

wherein the instructions for estimating the ray energy reflected by each sub-volume of 
the subset include 

estimating a maximum energy differential of the sub- volume; 

comparing the maximum energy differential against a predefined energy error 

threshold; 

if the maximum energy differential is above the predefined energy error 
threshold, recursively 

selecting a smaller sub-volume along the ray path; and 
estimating a new maximum energy differential of the smaller sub- 
volume; and 

if the maximum energy differential is below the predefined energy error 
threshold, calculating the amount of ray energy reflected by the sub-volume using the 
illumination model. 
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34. (Original) The system of claim 33, wherein the maximum energy differential 
depends on the opacity transfer function and the maximum, average, and minimum data 
values of the sub-volume. 

35. (Original) The system of claim 33, wherein the amount of ray energy reflected 
by the sub-volume depends on the length of ray path within the sub- volume, the opacity 
transfer function within the sub-volume, the average scalar field value of the sub-volume, and 
the local gradient vector of scalar field within the sub-volume. 

36. (Original) The system of claim 33, wherein if the sub-volume is a smallest 
sub-volume comprising 2x2x2 3-D cells, the smaller sub-volume is a 3-D cell within the 
smallest sub-volume, and if the opacity function does not vary monotonically within the cell, 
the 3-D cell is further divided into multiple sub-cells until the dimension of a smallest sub- 
cell reaches the cross-section of the ray. 

37. (Original) The system of claim 33, wherein if the sub-volume is a smallest 
sub-volume comprising 2x2x2 3-D cells, the smaller sub-volume is a 3-D cell within the 
smallest sub-volume, and if the opacity function varies monotonically within the cell, the 
maximum energy differential of the 3-D cell is estimated by dividing the maximum energy 
differential of the sub-volume by 2. 

38. (Original) The system of claim 33, wherein if the sub-volume is a smallest 
sub-volume comprising 2x2x2 3-D cells, the smaller sub-volume is a 3-D cell within the 
smallest sub- volume, and if the opacity transfer function varies monotonically within the cell 
and an iso-surface exists in the 3-D cell, the maximum energy differential of the 3-D cell is 
calculated using the eight data values at the corners of the 3-D cell and the opacity transfer 
function. 

39. (Original) The system of claim 33, wherein the predefined energy error 
threshold is modulated by an image rendering speed specified by a user, and a zoom factor in 
the case of parallel projection or a perspective angle and a perspective distance between the 
image plane and the 3-D dataset in the case of perspective projection. 

40. (Original) The system of claim 2 1 , wherein 

the instructions for estimating pixel values at other locations on the 2-D image 

plane include 

for each location, 
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selecting four pixel values associated with four ray origins surrounding 

the location; and 

bi-linearly interpolating a pixel value at the location using the four 

pixel values. 

41 . (Currently Amended) A computer program product for use in conjunction with 
a computer system, the computer program product comprising a computer readable storage 
medium and a computer program mechanism embedded therein, the computer program 
mechanism comprising: 

instructions for fragmenting a sampled 3-D dataset of a scalar field into a 
plurality of sub-volumes of different sizes, each sub- volume associated with a set of data 
value parameters characterizing the data value distribution of the scalar field within the sub- 
volume; 

instructions for defining an opacity transfer function that is dependent upon 
data values of the scalar field and an illumination model; 

instructions for selectively casting a plurality of rays from a 2-D image plane 
towards the sampled dataset, each ray having an initial ray energy and a cross-section; 

for each ray launched from a selected location on the 2-D image plane, 

instructions for selecting a subset of the plurality of sub- volumes for 
interacting with the ray; 

instructions for estimating the ray energy reflected by each sub-volume 
of the subset using the opacity transfer function and the illumination model; and 

instructions for summing the reflected ray energy as a pixel value at 
the selected location on the 2-D image plane; and 

instructions for estimating pixel values at other locations on the 2-D image plane 
using the pixel values at the selected locations 

wherein the instructions for selectively casting a plurality of ravs from a 2-D 
image plane include 

subdividing the 2-D image plane into a plurality of sub-planes; 
for each of the plurality of sub-planes, 

casting four ravs from the four corners of the sub-plane and estimating 
a pixel value at each corner; 

calculating a maximum pixel value variation within the sub-plane; and 
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recursively subdividing the sub-plane into multiple child sub-planes of 
smaller sizes bv casting a ray from the center of the sub-plane until the maximum pixel value 
variation of the sub-plane is below a predefined imaging error threshold . 

42. (Original) The computer program product of claim 41, wherein 
the step of fragmenting a sampled 3-D dataset includes 

fragmenting the 3-D dataset into eight sub-volumes; and 
for each sub-volume, recursively fragmenting it into eight smaller sub- 
volumes until the size of the smallest sub-volumes reaches a predefined size limit. 

43. (Original) The computer program product of claim 42, wherein the predefined 
size limit is a sub-volume comprising 2x2x2 3-D cells and the eight corners of each cell are 
associated with eight data values of the scalar field. 

44. (Original) The computer program product of claim 43, wherein the data value 
at any location within the cell is tri-linearly interpolated using the eight data values at the 
eight corners of the cell. 

45. (Original) The computer program product of claim 41 , wherein the set of 
parameters include a maximum, an average value, and a minimum data value of the scalar 
field within the sub- volume. 

46. (Original) The computer program product of claim 41 , further comprises 
instructions for constructing an octree comprising a root node, a plurality of 

intermediate nodes, and a plurality of leaf nodes; 

instructions for associating the root node with the 3-D dataset; 

instructions for associating each of the plurality of leaf nodes with a smallest sub- 
volume from the plurality of sub-volumes; and 

instructions for associating each of the plurality of intermediate nodes with a sub- 
volume from the plurality of sub- volumes that is larger than the smallest sub- volume. 

47. (Canceled) 

48. (Currently amended) The computer program product of claim [[47]] 41, 
wherein the maximum pixel value variation within the sub-plane is defined as the maximum 
deviation of pixel values at the four corners of the sub-plane from an average pixel value of 
the sub-plane. 

49. (Currently amended) The computer program product of claim [[47]] 41, 
wherein the predefined imaging error threshold is modulated by an image rendering speed 
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provided by a user, a distance to an edge of an object embedded in the 3-D dataset, and a 
difference between a pixel value estimated from an adaptive ray casting and a pixel value 
estimated from a bi-linear interpolation. 

50. (Currently amended) Th e computer program product of claim 11, A computer 
program product for use in conjunction with a computer system, the computer program 
product comprising a computer readable storage medium and a computer program 
mechanism embedded therein, the computer program mechanism comprising: 

instructions for fragmenting a sampled 3-D dataset of a scalar field into a 
plurality of sub-volumes of different sizes, each sub-volume associated with a set of data 
value parameters characterizing the data value distribution of the scalar field within the sub- 
volume; 

instructions for defining an opacity transfer function that is dependent upon 
data values of the scalar field and an illumination model; 

instructions for selectively casting a plurality of rays from a 2-D image plane 
towards the sampled dataset, each ray having an initial ray energy and a cross-section; 

for each ray launched from a selected location on the 2-D image plane, 

instructions for selecting a subset of the plurality of sub-volumes for 
interacting with the ray; 

instructions for estimating the ray energy reflected by each sub-volume 
of the subset using the opacity transfer function and the illumination model; and 

instructions for summing the reflected ray energy as a pixel value at 
the selected location on the 2-D image plane; and 

instructions for estimating pixel values at other locations on the 2-D image 
plane using the pixel values at the selected locations 

wherein the instructions for selecting a subset of the plurality of sub-volumes 
for interacting with the ray include 

identifying a largest sub-volume along the ray path and its corresponding 
maximum and minimum data values; 

checking if the opacity transfer function varies monotonically between the 
maximum and minimum scalar field values; 

if the function does not vary monotonically, recursively 
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identifying a smaller sub-volume along the ray path and its 
corresponding maximum and minimum data values; and 

checking if the opacity transfer function varies monotonically between 
the maximum and minimum scalar field values of the smaller sub-volume; and 

if the function does vary monotonically, calculating the amount of ray energy 
reflected by the sub-volume during its interaction with the ray. 

5 1 . (Original) The computer program product of claim 50, wherein two lookup 
tables are constructed for the opacity transfer function such that a forward lookup table 
contains the data value difference to a nearest local extreme of the opacity transfer function 
along the data value increasing direction and a backward lookup table contains the data value 
difference to a nearest local extreme of the opacity transfer function along the data value 
decreasing direction. 

52. (Original) The computer program product of claim 51, wherein if the 
maximum data value of the sub-volume is smaller than the summation of the minimum data 
value of the sub-volume and its corresponding data value difference stored in the forward 
lookup table or the minimum data value of the sub- volume is larger than the difference 
between the maximum data value of the sub-volume and its corresponding data value 
difference stored in the backward lookup table, the opacity transfer function varies 
monotonically between the minimum and maximum data values. 

53 . (Currently amended) Th e comput e r program product of claim M, A computer 
program product for use in conjunction with a computer system, the computer program 
product comprising a computer readable storage medium and a computer program 
mechanism embedded therein, the computer program mechanism comprising: 

instructions for fragmenting a sampled 3-D dataset of a scalar field into a 
plurality of sub-volumes of different sizes, each sub-volume associated with a set of data 
value parameters characterizing the data value distribution of the scalar field within the sub- 
volume; 

instructions for defining an opacity transfer function that is dependent upon 
data values of the scalar field and an illumination model; 

instructions for selectively casting a plurality of rays from a 2-D image plane 
towards the sampled dataset, each ray having an initial ray energy and a cross-section; 

for each ray launched from a selected location on the 2-D image plane, 
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instructions for selecting a subset of the plurality of sub- volumes for 
interacting with the ray; 

instructions for estimating the ray energy reflected by each sub-volume 
of the subset using the opacity transfer function and the illumination model; and 

instructions for summing the reflected ray energy as a pixel value at 
the selected location on the 2-D image plane; and 

instructions for estimating pixel values at other locations on the 2-D image 
plane using the pixel values at the selected locations 

wherein the instructions for estimating the ray energy reflected by each sub- 
volume of the subset include 

estimating a maximum energy differential of the sub-volume; 

comparing the maximum energy differential against a predefined energy error 

threshold; 

if the maximum energy differential is above the predefined energy error 
threshold, recursively 

selecting a smaller sub-volume along the ray path; and 
estimating a new maximum energy differential of the smaller sub- 
volume; and 

if the maximum energy differential is below the predefined energy error 
threshold, calculating the amount of ray energy reflected by the sub-volume using the 
illumination model. 

54. (Original) The computer program product of claim 53, wherein the maximum 
energy differential depends on the opacity transfer function and the maximum, average, and 
minimum data values of the sub-volume. 

55. (Original) The computer program product of claim 53, wherein the amount of 
ray energy reflected by the sub-volume depends on the length of ray path within the sub- 
volume, the opacity transfer function within the sub-volume, the average scalar field value of 
the sub-volume, and the local gradient vector of scalar field within the sub-volume. 

56. (Original) The computer program product of claim 53, wherein if the sub- 
volume is a smallest sub-volume comprising 2x2x2 3-D cells, the smaller sub-volume is a 3- 
D cell within the smallest sub-volume, and if the opacity function does not vary 
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monotonically within the cell, the 3-D cell is further divided into multiple sub-cells until the 
dimension of a smallest sub-cell reaches the cross-section of the ray. 

57. (Original) The computer program product of claim 53, wherein if the sub- 
volume is a smallest sub-volume comprising 2x2x2 3-D cells, the smaller sub-volume is a 3- 
D cell within the smallest sub- volume, and if the opacity function varies monotonically 
within the cell, the maximum energy differential of the 3-D cell is estimated by dividing the 
maximum energy differential of the sub- volume by 2. 

58. (Original) The computer program product of claim 53, wherein if the sub- 
volume is a smallest sub-volume comprising 2x2x2 3-D cells, the smaller sub-volume is a 3- 
D cell within the smallest sub-volume, and if the opacity transfer function varies 
monotonically within the cell and an iso-surface exists in the 3-D cell, the maximum energy 
differential of the 3-D cell is calculated using the eight data values at the corners of the 3-D 
cell and the opacity transfer function. 

59. (Original) The computer program product of claim 53, wherein the predefined 
energy error threshold is modulated by an image rendering speed specified by a user, and a 
zoom factor in the case of parallel projection or a perspective angle and a perspective 
distance between the image plane and the 3-D dataset in the case of perspective projection. 

60. (Original) The computer program product of claim 4 1 , wherein 

the instructions for estimating pixel values at other locations on the 2-D image plane 

include 

for each location, 

selecting four pixel values associated with four ray origins surrounding the 

location; and 

bi-linearly interpolating a pixel value at the location using the four pixel 

values. 

61.-68. (Canceled) 



1 -NY/1 998980.1 



- 18- 



